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ABSTRACT 

Context. During the last 30 years the surface of Pluto has been characterized, and its variability has been monitored, through contin¬ 
uous near-infrared spectroscopic observations. But in the visible range only few data are available. 

Aims. The aim of this work is to define the Pluto’s relative reflectance in the visible range to characterize the different components of 
its surface, and to provide ground based observations in support of the New Horizons mission. 

Methods. We observed Pluto on six nights between May and luly 2014, with the imager/spectrograph AC AM at the William Herschel 
Telescope (La Palma, Spain). The six spectra obtained cover a whole rotation of Pluto (Proi = 6.4 days). For all the spectra we 
computed the spectral slope and the depth of the absorption bands of methane ice between 0.62 and 0.90 ^m. To search for shifts of 
the center of the methane bands, associated with dilution of CH 4 in N 2 , we compared the bands with reflectances of pure methane ice. 
Results. All the new spectra show the methane ice absorption bands between 0.62 and 0.90 pm. The computation of the depth of the 
band at 0.62 pm in the new spectra of Pluto, and in the spectra of Makemake and Eris from the literature, allowed us to estimate the 
Lambert coefficient at this wavelength, at a temperature of 30 K and 40 K, never measured before. All the detected bands are blue 
shifted, with minimum shifts in correspondence with the regions where the abundance of methane is higher. This could be indicative 
of a dilution of CH4:N2 more saturated in CH4. The longitudinal and secular variations of the parameters measured in the spectra are 
in accordance with results previously reported in the literature and with the distribution of the dark and bright material that show the 
Pluto’s albedo maps from New Horizons. 

Key words. Kuiper belt objects: individual: Pluto - methods: observational - methods: numerical - techniques: spectroscopic - 
planets and satellites: composition 


1. Introduction 


The surface of Pluto is covered with ices of varying degrees of 
volatility, and even at its low average temperature of ~40 K, there 
is an exchange of molecules between the surface and the thin at¬ 
mosphere. Because of Pluto’s high obliquity (119.6°) and the ec¬ 
centricity of its orbit (perihelion 29.66 AU, aphelion 48.87 AU), 
there are pronounced seasons over the 247.7-year orbital period. 
These factors ensure that Pluto’s surface undergoes changes on 
a seasonal time scale. Secular changes may also occur, owing 
to the long-term evolution of the chemical components of the 
surface and atmosphere under the influence of solar and cosmic 
ray irradiation. On the time scale of Pluto’s 6.4-day rotation pe¬ 
riod, clear changes in the distribution of ices across its surface 
are seen. 


Indeed, on seasonal and rotational time scales, changes in 
Pluto’s overall (global) brightness and in the near-infrared spec¬ 
tral signatures of the three princ ipal ices, N2, CH4 , and CO 
have b een detected and quanti fied dOrundv et al.ll20ll (hereafter 
G13); ICruikshank et al.ll2015l and references therein). Because 
Pluto has only been observed over one-third of its full orbit, sec¬ 
ular changes might be confused with seasonal variations, with 


the result that at the present time there is no unequivocal evi¬ 
dence for observable secular changes. 

In addition to the ices on Pluto’ surface, there is another ma¬ 
terial present that imparts colors as seen in the extended visible 
spectral region, 0.30-1.00 pm. As deduced from multi-color pho¬ 
tometric measurements and images, th e colors vary from light 
to medium yellow, red-brown, to black (iBuie et alJl2010bl) . The 
ices detected to date do not produce these colors, and it is gener¬ 
ally accepted that photochemical reactions in Pluto’s atmosphere 
and in the surface ices create colored macromolecular matter that 
contains or is made of complex orga nic chemicals derived from 
carbon, nitrogen, and oxygen (e.g.. iMaterese et ani2014l) . The 
coloring properties of these materials are well established (e.g. , 
Khare_etj^ 1^84^ Cmikshank et alj|2005t llmanaka et al.ll201^ 


Materese et al.ll201^ . The coloring agent is not uniform over 

Pluto’s surface, changing noticeably over the dwarf-planet’s 6.4- 
day rotation. 

The best way to quantify the color of Pluto and to search 
for diagnostic characteristics is to record the reflectance proper¬ 
ties of the surface in the photo-visual, or extended visual spec¬ 
tral region (0.30 - 1.00 pm). Only a few such measurements 
have been reported in the literature, beginning with observations 
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made in 1970 by Fixe]^^ (Il970l) . and followed bv iBenner et al.l 
(fT978h . lBellet al.1 (1197^ fdata pu bl ished in Cruikshank and 
Brown 1986), Barker et all (|1980l) . iBuie & FinkI (Il987h and 
iGrundv & FinkI (Il996l) . Individually or together, these studies 
do not dehne a well sampled range of rotation and seasonal 
change, or a full sample of wavelengths in the desired range at 
adequate spectral resolution. The history of color measurements 
with ground-based tel escopes and the Hubble S pace Telescope 
(HST) is reviewed by ICruikshank et al.l (12015h . who note that 
the color measured in the (B-V) photometric index was constant 
from 1953 to 2000 withi n a few tenths of a percent, but then 
changed abruptly in 2002. [Buie et alJ (l2010bh give the details of 
the HST observations that showed a global change in color to a 
redder (B-V) index that persists to the present day. 

The variation of the relativ e amounts of ices a nd dark mate¬ 
rial have been mapped before (iBuie et al.ll2010bt G13), but are 
even more obvious now in the light of the preliminary new im¬ 
ages obtained by the New Horizons mission. The maps combin¬ 
ing the reflectance and colours of Pluto’s surface clearly show a 
diverse palette from the very dark region on the equator, infor¬ 
mally named Cthulhu Regio to the “clear peach orange” region, 
informally named the Tombaugh Regio (Stern et al., submitted). 

During 2014, we conducted a new observational study in or¬ 
der to provide a high-quality data set of spectrophotometric data 
in the 0.40 - 0.93 pm spectral range. These observations en¬ 
able the dehnition of Pluto’s reflectance and the investigation 
of weak methane ice absorption bands shortward of 1.00 pm. 
The data have hemispheric spatial resolution over a whole rota¬ 
tion of the Pluto’s longitudes. The new data are compared with 
earlier observations to detect any changes in spectral reflectance 
that may have occurred in the last few decades. This comparison 
also provides a basis for the modeling of the spectrum with ra¬ 
diative transfer models that account for the known components 
of Pluto’s surface as well as additional components that are sus¬ 
pected to occur. The observations presented here include the re¬ 
flected light of both Pluto and Charon. These observations are 
also in support of the New Horizons mission that observed the 
Pluto system in July, 2015 and are part of a long term campaign 
to follow the secular evolution of the Pluto’s surface properties. 

In the next section we describe the observations. In section 3 
we analyze the obtained spectra, calculating the slopes of the rel¬ 
ative reflectances, the depth and the area of the absorption bands 
of methane ice, and measuring the shifts of their centers. In Sec¬ 
tion 4 we discuss the results and in Section 5 we report our con¬ 
clusions. 


2. The observations 

We observed Pluto on six different nights between May 13 and 
July 17, 2014 using the low-re solution spectros copic mode of 
ACAM (auxiliary-port camera, feenn et al]l2008h . ACAM is an 
imager/spectrograph that is mounted permanently at a folded- 
Cassegrain focus of the 4.2 m William Herschel Telescope at the 
Observatorio del Roque de los Muchachos (ORM), La Palma. 
We used the V400 disperser and the 1.5” slit providing a resolu¬ 
tion of ~300 at 0.55 pm (the dispersion is 3.3 A/px). The slit was 
oriented at the parallactic angle, and the tracking of the telescope 
was set to compensate for Pluto’s proper motion. 

In order to cover a whole rotation of Pluto, we originally 
scheduled seven observations equally spaced in time to observe 
Pluto’s surface every ~50 degrees in longitude. Unfortunately, 
we discarded one of those observations due to technical prob¬ 
lems. The six that we kept cover 310 degrees in longitude. The 
details of the observations are described in Table [T] We define 


longitudes following the IAU2009 model that corresponds with 
the right-hand rule. Every night we obtained three consecutive 
spectra of Pluto of 300 seconds exposure time each, wich pro¬ 
vides a S/N > 200 in most of the wavelength range and even 
better after averaging the three spectra. 

Data reduction was performed using the Ima ge Reductio n 
and Analysis Facility (IRAF) standard procedures (lTodvlll993h . 
Preprocessing of the CCD images included bias and flat field cor¬ 
rection. Flat field images were obtained using calibration lamps 
(tungsten continuum lamp) in the afternoon. Wavelength cali¬ 
bration lamps (CuAr and CuNe) were obtained before and after 
the set of Pluto’s spectra, maintaining the same position angle to 
avoid the effect of instrumental flexures in the wavelength cali¬ 
bration. 


To correct the spectra for telluric absorptions and to obtain 
the relative re flectance, we observed every night two different 
Landolt stars ( Landolt 199^ that have proven to be good solar 
analogs (iLicandro et alJ 20^1) . We observed them at airmasses 
very similar to that of Pluto (see Table [T] for the details). These 
solar analogs are bright (10 < Vmag < 12.4) so the exposure 
times were always shorter than 60 seconds. To reduce the num¬ 
ber of bad pixels and artifacts, we obtained three spectra of each 
of the stars as we did for Pluto. The reduction of the spectra of 
the solar analogs was done following the same procedure as for 
Pluto. All the calibration images, the solar analogs and the Pluto 
images were obtained in ~1 hour. 

The extinction of the sky is dependent on the wavelength, 
with shorter wavelengths experience greater extinction. To mini¬ 
mize the spectral extinction effect from the difference in airmass 
between the stars and the target, we applied color correction to 
the spectra of the object and the stars. We used the mean extinc- 
tion coeffici ent for the Observatorio del Roque de los Muchachos 
(lKinglll98.5h . 


We Anally divided the color-corrected spectra of Pluto by 
the color-corrected spectra of the stars to obtain Pluto’s relative 
reflectance. We used the sky lines at this step to refine the wave¬ 
length calibration. Before averaging all of them into a ‘flnal’ rel¬ 
ative reflectance, we visually inspected the 18 reflectances and 
discarded those that showed problems due to the systematic er¬ 
rors. We averaged the remainder of them into a flnal relative re¬ 
flectance for each date of observation. 


In Fig. [T] we present the six reflectance spectra obtained in 
this observing campaign. These are the average of the Pluto -i- 
Charon spectra divided by the spectra of the solar-type stars. This 
division cancels the solar color and the solar and telluric absorp¬ 
tion lines. The spectra have been normalized to 1.0 at 0.75 pm, to 
enhance the differences in the slope of the continuum below 0.70 
pm. The contribution of Charon has not been removed in this 
plot. Nevertheless, the spurious light of Charon is not a problem 
as Charon is two magnitudes fainter than Pluto, and its color is 
neutral at these wavelengths (as we explain in detail later). 

In Fig. |2] we present the data from Fig. [T] smoothed using 
an average Alter of 5 pixels and plotted in geometric albedo. To 
transform from reflectance to al bedo, we u s ed the magnitudes 
of Pluto alone from Fig. [T] of iBuie et al.1 (1201 Oal) . in which 
they plot the lightcurve of Pluto resolved from Charon. The 
values in the figure are shown as mean opposition magnitudes 
correcte d to phase angle a - 1°. We have taken the iBuie et al.l 
(1201 Oah magnitudes at a = 1° for the longitudes we observed 
and have calculated the geometric albedos corresponding to 
those longitudes. Using a radius of 1184 km for the dwarf 
planet (iLellouch et al.lUOTsI) . we calculated the global geometric 
albedo of Pluto alone at the V wavelength (0.56 pm) for each 
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Table 1. Details of the observations. Time = mid exposure time; am(Pluto) and am(SA) = mid exposure airmass of Pluto and the solar analogs, 
respectively; Lon = Pluto sub-Earth longitude; Lat = Pluto sub-Earth latitude; p,. = Geometric albedo of Pluto at 0.56 /rm. 


2014 Date 

Time 

am (Pluto) 

Solar Analog 

am (SA) 

Lon (deg) 

Lat (deg) 

P> 

May 13 

04:38 

1.5 

L112-1333 

1.6 

103.7° 

51.2° 

0.44 




L107-998 

2.3 




May 14 

03:40 

1.6 

L112-1333 

1.7 

49.7° 

51.2° 

0.48 




L107-998 

1.5 




June 6 

01:50 

1.7 

L112-1333 

1.5 

197.9° 

50.9° 

0.59 




L107-998 

1.2 




June 7 

01:57 

1.6 

L112-1333 

1.5 

141.3° 

50.9° 

0.50 




L107-998 

1.3 




June 18 

00:23 

1.8 

L112-1333 

2.3 

245.2° 

50.7° 

0.57 




P041C 

1.5 




July 17 

22:40 

1.7 

L112-1333 

1.7 

358.9° 

50.0° 

0.52 




P041C 

1.4 






Fig. 1. Relative reflectances of Pluto corresponding to hemispheres cen¬ 
tered at six different sub-Earth longitudes. All the reflectances are nor¬ 
malized at 0.60 yum and shifted by 0.15 in relative reflectance for clar¬ 
ity. The labels are ordered (top to bottom) in the same order than the 
reflectances. The vertical dashed lines show the telluric absorptions. 
Some of these telluric absorptions have not been completely removed. 
The absorption at 0.52 pm is an artifact, it is present on the spectra of the 
solar analogs, its center corresponds to the magnesium triplet (0.5169 
-0. 5183 pm) that is present in G stars. 


date we observed (listed in Table [TJ using the equation 

-2.5 log Pv = Ysun - yPluto -5 log r H- 5 log (RxA), 

where VjMn = -26.75, r = Pluto’s radius (in AU), R = Pluto’s 
mean heliocentric distance (39.482 AU) and A= Pluto’s geo¬ 
centric distance (38.482 AU). A small revision of Pluto’s radius 
from the New Horizons investigation makes a negligible change 
to the albedo calculations. 

We made no correction for the contribution of Charon to 
the color of (Pluto + Charon) that we observed, assuming that 
Charon’s color is neutral gray. From observation s of the mutua l 
transit and occultation events of the mid-1980s. iBinzej (Il988h 
found a nearly neutral gray color for Charon , with B-V = 0.700 
+ 0.010 ((B-V)sMn = 0.65). More recently. iBuie et aO (l2010ah 



Wavelength {fim) 

Fig. 2. Same spectra of Fig. [T] smoothed using an average filter of 5 
pixels and represented in geometric albedo. The symbol plus indicates, 
for each spectrum, the global geometric albedo of Pluto alone at 0.56 
pm calculated for the corresponding date of the observation. 

derived (B-V) = 0.7315 + 0.0013 using HST dat a. At the pho¬ 
tometric B band (0.44 uml. lTholen & Bul^ (Il990l) give the geo¬ 
metric albedo of Pluto as 0.44-0.61, and that for Charon 0.38. 

3. Results and analysis 

For each spectra we computed four parameters: (1) the spectral 
slope in order to quantify the color of Pluto at different lon¬ 
gitudes and, thus, search for possible variations in the amount 
of solid complex organics over its surface; ( 2 ) the integrated 
area of the absorption bands of methane ice at 0.62 pm in 
order to determine if this weak CH 4 band is detected as in 
other trans-Neptunian obje cts (TNOs) (iLicandro et al.l l^06bt 
lAlvarez-Candal et al.ll201 ih : (3) the depth of the methane ice 
bands between 0.62 and 0.90 pm in order to search for possible 
variations of the abundance of CH 4 -ice on the surface of Pluto 
with longitude; (4) the shift of the center of the absorption bands 
of methane ice at 0.73, 0.79, 0.80, 0.84, 0.87 and 0.89 pm in or¬ 
der to search for possible variations o f the dilution of CH 4 in N ? 
with longitude (previously reported bv lOuirico & Schmitl)l997l) . 
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Table 2. Estimated systematic error in the slope expressed in %/1000 
A. 


2014 Date 

Syst. error 

2014 Date 

Syst. error 

May 13 

0.45 

June 7 

0.76 

May 14 

0.70 

June 18 

0.73 

June 6 

0.37 

July 17 

0.71 


3.1. Slope 


The slope in the visible is an indicator of the amount, particle 
size, or type of complex organics on the surface. Complex organ¬ 
ics have been det ected on the surface of asteroids (24) Themis 
and (65) Cybele dRivkin & Emervl[201(A iLicandro et al.ll2Ql ih 
and the satu rnian moons lapetus and Phoebe ( Cruikshank et al.l 
l2008l |2014|) . and their presence has been i nferred for other 
solar system objects (e.g. TNOs and Titan, ICruikshank et alJ 
I2OO5I) . The recent discovery of a suite of 16 organic compounds, 
comprising numerous carbon and nitrogen-rich molecules, 
on the surface of the comet 67P/Chu ryomov-Gerasimenko 
dWright et al.ll2015llGoesmann et al.l2015h suggests that organic 
compounds that are the precursor of the more complex organics 
(such as tholins) are a ubiquitous ingredient of the minor bod¬ 
ies in the Solar System. Not surprisingly, red slopes, typically 
associated with the complex organics, are common among the 
spectra of TNOs. 

Here we compute the slope of Pluto’s visible spectrum in two 
ranges: from 0.40 to 0.70 pm and from 0.57 to 0.70 pm. The first 
range comprises the full spectral region below 0.70 pm covered 
by our spectra. The second range enables the comparison with 
previous spectra of Pluto, Eris, and Makemake in the literature. 
To compute the slope values in both cases, we first normalized 
all the spectra at 0.60 pm, and then we fitted each of them in the 
chosen range with a linear function. We made this computation 
for all the six spectra and for the average relative reflectance. The 
slope is expressed here as %/1000 A. 

To evaluate the errors in the computation of the slope, we 
took into account both the error in the fit and the systematic er¬ 
rors in the reflectance. The systematic errors depend on different 
contributions: e.g., errors introduced by color corrections or by 
dividing by the spectra of the solar analogs. The error introduced 
by dividing by the spectra of the solar analog is the dominant part 
of the systematic errors. To estimate this error we compared the 
six individual spectra of the s tars observed each n ight follow¬ 
ing the procedure described in lLorenzi et alJ (1201 4l) and Pinilla- 
Alonso (submitted). After discarding those spectra with an ex¬ 
treme behavior, we took the standard deviation of the slope of 
these spectra as the error in the slope for the night. These values 
are shown in Table |2l We take the maximum value in the run 
as the systematic error. We found that the formal error coming 
directly from the computation of the slope is an order of magni¬ 
tude lower than the systematic error, so we took 0.8 %/1000 A 
as the final uncertainty associated to the slope calculation for all 
the spectra. 

To compare with other spectra in the literature, we computed 
the slope of those spectra in the same way that we did for the new 
data presen ted herein. We u s ed for this comparison the spectra of 
Pluto from iGrundv & FinkI (Il996h and the spectra of two other 
large TNOs, (136472) Makemake from Licandro et al. (2006a) 
and (136199) Eris from lAlvarez-Candal et alJ (1201 Ih . For the 
spectra from the literature of Pluto, Makemake, and Eris, there 
are no estimates for the systematic errors in the original papers, 
so we take the typical value of 1 %/1000 A (e.g. lLazzaro et al.) 



Wavelength (^m) 


Fig. 3. Averaged spectrum of the six spectrum presented in this work, 
normalized at 0.60 pm. The positions of the absorption bands of CH 4 
are represented by dotted lines. 


l2004h for all the spectra (for these spectra the systematic eiTors 
are also an order of magnitude greater than those associated with 
the fit calculation). The results are listed in Table|3] 

3.2. Methane absorption bands 

The new spectra clearly show methane ice absorption bands at 
0.73, 0.78 - 0.80, and 0.83 - 0.91 pm. In Fig. [3j we show the 
spectrum of Pluto obtained averaging the six spectra presented 
in this work, with the methane ice absorption bands evidenced 
as vertical dotted lines. All these bands have very low absorp¬ 
tion coefficients, 5.0x10“^, ~1.3-1.4x10“^, ~ 1.0-8. Ox 10~^ cm~\ 
respe ctively, as measured in the laboratory at 40 K (iGrundv et al.l 
l2002h . This average spectrum also shows an absorption around 
0.62 pm. This band has been also observed in the spectra of other 
large TNOs dominated by methane ice, like Eris and Makemake 
(lAlyarez-Candal et^l2()l U ILicandro et alj|200^ see Fig. lU, 
so its identification is not in question. Howeyer, the absorption 
coefficient of this band is unknown as it has not been measured 
in the laboratory. If we look at each particular reflectance ob¬ 
tained in 2014, some of them show it yery clearly, in particularly 
the phase centered at 358.9°. 

To confirm our detection, we compute the integrated area of 
the absorption band. For this computation, we first fit a straight 
line to the continuum on both sides of the band, between 0.593 
- 0.615 pm and 0.624 - 0.653 pm, and diyide the spectrum by 
this continuum (see Fig. |5]l. A first look at the plots in Fig. |5] 
suggests a detection of the band, although in some of the phases 
(e.g. 49.7° or 103.7°) it is not clear if the absorption is significant 
aboye the noise. 

We then calculate the area of the absorption band after diyi- 
sion by the continuum, integrating it between 0.615 and 0.624 
pm. We propagate the errors, taking into account in each point 
the error associated with the relatiye reflectance and the error 
in the calculation of the fit of the continuum. Finally, we iter¬ 
ate the computation of the area, changing randomly the limits of 
the band, and we take the mean yalue and the standard deyia- 
tion as the final yalues for the area and the associated error. No¬ 
tice that for all the spectra (eyen those where the detection was 
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Table 3. Computed spectr al slope between 0 . 57 and 0.70 fxm, and between 0.40 and 0.70 /rm, with the spectra normalize d at 0.60 fim, for the 
spectra of Pluto (this work ;|Buie & Finkll98^ : lGrundv & FinkllQQhh . for the spectrum of Makemake dLicandro et al.l2006bh and for the spectrum 
of Eris dAlvarez-Candal et alJl2011h . 


Longitude (degrees) 

Slope (%/1000 A ) 0.57-0.70 pm 

Slope (%/1000 A ) 0.40-0.70 pm 

Reference 

49.7° 

14.2 

19.9 

This work 

103.7° 

12.6 

20.6 


141.3° 

12.5 

19.3 


197.9° 

12.7 

19.4 


245.2° 

10.7 

17.3 


358.9° 

10.8 

17.2 


Average spectrum 

12.3 

19.0 


25°(1983) 

5.7 


Buie & Fink (19871 

329° (1983) 

5.1 

- 


144° (1990) 

7.7 

- 

Grundv & Fink (19961 

56° (1993) 

10.8 

- 

5? 

150° (1994) 

10.4 

- 


Makemake 

12.9 

17.3 

Licandro et al. (200661 

Eris 

4.4 

6.5 

Alvarez-Candal et al. 120111 



Wavelength (^m) 


Fig. 4. Comparison of the spectra of Eris, Pluto and Makemake. The 
three spectra are normalized at 0.60 ^m. All three objects show the ab¬ 
sorption bands of methane ice. The band at 0.62 ;um is more evident on 
Makemake and Eris than on Pluto. Note the similarity between the slope 
of Pluto and Makemake that suggests the presence on both surfaces of 
a very similar coloring agent. 


suspicious) the integrated area of the band is a positive value. 
Some of the integrated areas, however, have large uncertainties. 
To confirm the detection and discard a serendipitous result, we 
do a second verification. 

For this test we repeat the procedure on a region of the con¬ 
tinuum where no band is expected. We define this continuum in 
a region with the same number of points between 0.55 and 0.57 
pm. What we could expect from this test is that we obtain both 
positive and negative values for the area, and that the absolute 
values of the areas would be smaller than the uncertainty. We 
show in Table|4]the results of this test. Notice that for the calcu¬ 
lation on the continuum, we obtained both positive and negative 
values, whereas for the absorption bands at 0.62 pm we obtained 







0.60 0.61 0.62 0.63 0.64 
Wavelength (^m) 


0.60 0.61 0.62 0.63 0.64 
Wavelength {^m) 


Fig. 5. Continuum-removed in the region of the spectra around the 0.62 
pm absorption band of methane. In red the absorption bands, in blue the 
continuum. 
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Table 4. Integrated area in unit of of the absorption band at 0.62 
pm and test on the continuum between 0.55 and 0.57 pm measured as 
explained in the text. 


Longitude 

Band at 0.62 pm 

Cont. 0.55 - 0.57 pm 

49.7° 

0.31+0.30 

0.07 + 0.70 

103.7° 

0.39 ± 0.21 

0.51+0.54 

141.3° 

0.29 + 0.13 

0.30 + 0.35 

197.9° 

0.33 + 0.14 

-0.11 +0.41 

245.2° 

0.37 + 0.14 

0.26 + 0.46 

358.9° 

0.66 + 0.16 

-0.57 + 0.31 

average spectrum 

0.35 + 0.09 

0.13 + 0.17 


all positive values. After this test, we can consider the detection 
of the absorption bands at 0.62 pm as reliable, with the phase 
centered at ~360° as the most clear detection. 

3.3. Band depth 

To evaluate the globally averaged distribution of CH 4 -ice around 
Pluto’s surface, we computed the depth of the absorption bands 
as D = l-R^j/Rf, where R* is the reflectance in the center of the 
band, and R^ the reflectance of the continuum at the same wave¬ 
length. As several of the absorption bands of CH 4 in our spectra 
are a superposition of different absorption bands, it is not possi¬ 
ble to measure the integrated area or the equivalent width of each 
single band. We calculated the band depth that can be isolated; 
0.62, 0.73, 0.79, 0.80, 0.84, 0.87, 0.89 and 0.90 pm for our six 
spectra of Pluto. To evaluate the uncertainty in the calculation of 
the depth we propagated the errors, taking into account the un¬ 
certainties associated with the value of the relative reflectance in 
the center of the band and with the value of the continuum at the 
same wavelength. In Table |5] we list the values obtained for the 
depths (in %) of the six spectra of Pluto and the corresponding 
errors. 


3.4. Lambert Coefficients 


The appearance of an intrinsically weak band like the one at 0.62 
pm shown in Sec. 13.31 demonstrates a very long (several cm) 
path-length of sunlight in the nitrogen ice in which CH 4 is in¬ 
corporated as a trace constituent. This, in turn, helps to set the 
size scale for the individual elements of polycrystalline solid N 2 
characteristic of the surface of Pluto. 

A great effort has been made to estimate the absorption co- 
effici ents of methane ice by means of laboratory experiments 
(e.g. iGrundv et al.ll2002h . However, these experiments do not 
typically cover wavelengths below 0.70 pm, as a large slab of 
methane is needed to produce that absorption. We will use this 
band observed in the spectrum of Pluto, Eris, and Makemake to 
estimate the absorption coefficient of CH 4 ice at 0.62 pm. For 
this purpose, we plotted on a logarithmic scale the estimated 
band depth for the bands in the 0.73-0.90 pm range for the spec¬ 
tra of Pluto (for each band we used the average of the six values 
obtained for the six spectra) as a func tion of the correspon ding 
Lambert coefficients at T = 40 K from iGrundv et al.l (l2002h . Fi¬ 
nally we fit a 1 -degree polynomial (i.e. straight line) to the points 
with a least squares fit (see Fig.| 6 ]top panel). 

We do the same for the spectrum of Makemake from 
Licandro et al.l ( 2006bl) and for the spectrum of Eris from 
Alvarez-Candal et al.l fconT) . In the case of Eris, we used the 
Lambert coefficients corresponding to T = 30 K, as this tem¬ 





Fig. 6. Band depth (in natural logarithmic scale) vs Lambert coefficients 
(on natural logarithmic scale) for Pluto (top panel), Makemake (middle 
panel) and Eris (bottom panel). 


perature is more similar to that estimated for the surface of this 
object (see Fig. | 6 ]bottom panel). In both cases we measured the 
band depth of methane ice in the spectrum of Makemake and 
Eris from the literature following the same procedure described 
in the previous section for Pluto. 

In all the plots in Fig. | 6 ] the blue line is the least squares fit, 
the discontinued line corresponds to the measured band depth at 
0.62 pm, with its corresponding error (dotted lines). The inter¬ 
section between the line fitting the points and the measured band 
depth value for 0.62 pm band (red point) represents the extrapo¬ 
lated value of the Lambert coefficient for the 0.62 pm band, and 
its error (red line). In Table 6 we listed the band depths (%), with 
the respective errors, and the value of the Lambert coefficient for 
each band, taken from the literature. 

The value of the extrapolated Lambert coefficients at 0.62 
pm for Pluto is 2.88x10“^ + 1.23x10“^ cm“', while the value 
obtained for Makemake is 2.23x10“^ + 0.85x10“^ cm“*. The two 
values are in agreement within the errors as expected, since both 
objects have similar estimated surface temperature. 

The value obtained for Eris is 2.97x10^^ + 1.46x10“^ cm“'. 
The difference due to the temperature is not significant within 
the error in the calculation. 
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Table 5. Band depth (%) and one-sigma errors for eight methane bands at each of the six longitudes observed. 


Longitude 0.62 /rm 0.73 yum 0.79 yum 0.80 yum 0.84 yLtm 0.87 yum 0.89 yum 0.90 yum 


49.7° 

1.16 

+ 

0.96 

4.97 

+ 

0.74 

2.57 

+ 

1.23 

2.52 

+ 

1.10 

1.49 ± 

1.63 

7.07 

+ 

1.59 

17.35 + 

1.34 

10.19 + 0.81 

103.7° 

0.89 

+ 

0.74 

4.27 

± 

1.03 

0.97 

± 

0.58 

1 .12+ 

0.68 

1.54 + 

0.59 

5.69 

± 

1.64 

16.03 + 

1.40 

9.98+0.55 

141.3° 

0.50 

+ 

0.48 

4.08 

± 

0.90 

1.33 

± 

1.20 

1.50 

+ 

1.04 

1.59+ 

1.56 

5.13 

± 

1.98 

15.77+ 

2.34 

9.85 + 0.54 

197.9° 

0.75 

+ 

0.41 

5.16 

± 

0.46 

0.93 

+ 

0.57 

1.58 

+ 

0.53 

1.05+ 

0.46 

6.68 

± 

2.64 

16.74+ 

1.17 

9.89 + 1.44 

245.2° 

0.89 

+ 

0.67 

4.99 

± 

1.04 

1.77+ 

0.44 

2.37 

+ 

0.71 

1.90 + 

1.05 

6.28 

± 

1.87 

17.33+ 

1.53 

11.01+0.73 

358.9° 

1.13 

+ 

0.38 

6.29 

± 

0.86 

1.85 

± 

0.86 

3.00 

+ 

0.61 

1.75 + 

0.69 

8.47 

± 

1.49 

19.13 + 

1.41 

10.41 + 1.09 


Table 6. Band depths (%), with respective errors, and Lambert coefficients dGrundv et alJl200^ used for the extrapolation of the Lambert coef¬ 
ficient at 0.62 yum. The values of the coefficient in bold format have been estimated i n this work. * average of the values computed for the six 
spectra shown in Table [3 ** band depth measured for the spectrum of Makemake from iLicandro et ^ ll2006bh and for the spectrum of Eris from 
lAlvarez-Candal et all d^l Ih . 


Object 

0.62 yum 

0.73 yum 

0.79 jum 

0.80 yum 

0.84 jum 

0.87 yon 

0.89 yum 

0.90 yum 

Pluto * 

Makemake** 

Lambert Coef. (40 K) 

0.9 + 0.3 
10.0+ 1.3 

2.23-2.88x10 3 

5.0+ 0.8 
28.8 + 2.5 
5.03x10-2 

1.6+ 0.6 
17.1 + 3.8 
1.29x10-2 

2.0 + 0.7 
16.7 + 2.2 
1.41x10-2 

1.6+ 0.3 
19.1+5.8 
1.03x10-2 

6.6+ 1.2 
33.2 + 3.7 
7.06x10-2 

17.1+ 1.2 
58.3 + 7.3 
5.10x10-' 

10.2 +0.4 
39.4 + 4.7 
1.69x10-' 

Eris** 

Lambert Coef. (30 K) 

4.6 + 1.0 

2.97x10 3 

14.9+ 1.7 
5.22x10-2 

9.9+ 1.9 
1.42x10-2 

9.4+ 1.0 
1.51x10-2 

6.8+ 1.7 
1.12x10-2 

17.9 + 3.2 
7.45x10-2 

42.7 + 3.5 
5.50x10-' 

25.1 + 1.7 
1.79x10-' 


3.5. Shifts of the center of the bands 


4. Discussion 


The centers of the absorption bands of methane ice provide 
information about the m ixing ratio of methane and nitrogen 
dOuirico & Schmit3ll997h . When methane is diluted in nitro¬ 
gen ice in small proportions, the cen ters of the methane band s 
are shifted towards bluer wavelengths dOuirico & Schmitt! 1997h . 
and the shift is larg er as the cont ent of nitrogen increases 


(Brunette et al. I 2 OO 8 ). Some works ( 

Ouirico & Schmitt 1 19971 

lOuirico et al.lll999t Doute et al.l 11999 

; Merlin et al.l 

201 Oh dis- 


cuss the existence of regions formed of pure CH 4 segregated 
from the N 7 :CH 4 so l ution. This effect was first reported for 
Pluto by lOwen et^ dl993h . However, this is in conflict with 
both altern ative models fo r explaining Pluto’s atmosphere. Ac¬ 
cording to iTraftonI d2Q15h . and references therein, the surface 
of Pluto would instead be covered by a solid solution of these 
ices, N 2 :CH 4 , with different mixing ratios. Moreover, two dilu¬ 
tion regimes, one saturated in CH 4 and the other saturated in N 2 
would coexist on the surface of Pluto. The dilution regime satu¬ 
rated in CH 4 would be the responsible of the relatively unshifted 
part of the Pluto’s spectra, previously attributed to pure CH 4 . 


We measured the shifts of the center of the bands by com¬ 
paring our spectra with theoretical relative reflectances of pure 
methane ice. The syntheti c relative reflectances were calculated 
using a Shkuratov model (IShkuratov et al.lll999h and the opti cal 
constants of methane for T = 40 K from Grundy et al.l (l2002h . 


For comparison, we obtained a suit of relative reflectances 
considering different grain sizes, then we shifted each model 
in wavelength until we found the best comparison. We used a 
chi-square test to choose the best combination of grain size and 
shift. To evaluate the error associated with the measurement, we 
repeated the process several times, excluding randomly some 
points of our spectra and changing the limits of the bands. The 
obtained shifts are listed in Table [T] In some cases, due to the 
low signal to noise or to the presence of some artifacts that ap¬ 
pear during the reduction process and change the shape of the 
band, it has been impossible to measure the band shift. 


Previous works have shown that Plut o presents a surface 
with a zonal (longit udinal) heterogeneity (iGrundv & Finkl[r996l: 
iBuie et aDl2010allfl G13). This heterogeneit y is clearly seen in 
processed images of Pluto from the HST dBuie et al.l 1201 Obi) 
and in the recently released images of the New Horizons mis¬ 
sion (www.nasa.goy/mission_pages/newhorizons). The hetero¬ 
geneity is associated with an uneyen distribution of the ices 
(CH 4 , N 2 , CO) and the complex organics on the surface. 

G13 report long time trends in near-infrared spectra obtained 
between 2000 and 2012 from the analysis of the absorption 
bands of CH 4 , CO, and N 2 present in the spectra. They find that 
the absorptions of CO and N 2 are concentrated on Pluto’s anti- 
Charon hemisphere (~180°), unlike absorptions of less yolatile 
CH 4 ice that are offset by roughly 90 degrees from the longitude 
of maximum CO and N 2 absorption. In the time scale of their 
obseryations, they find yariations in the depth and in the ampli¬ 
tude of the diurnal yariation of the absorption bands of CH 4 , CO 
and N 2 . They conclude that the geometric distribution of these 
ices (in particular CO and N 2 ) on the surface of Pluto can ex¬ 
plain these secular changes, although seasonal yolatile transport 
could be at least partly responsible. 

Although yisible wayelength spectroscopy does not enable 
the study of the N 2 and CO concentrations, we can discuss the 
longitudinal yariation within our set of obseryations and the sec¬ 
ular yariation when comparing with preyious obseryations, of the 
distribution of methane-ice and of the dark-red organic material. 

Also, from the study of the center of the CH 4 bands, we can 
extract information on the relatiye amount of CH 4 diluted in N 2 . 

With respect to the methane, G13 find that the stronger CH 4 
absorption bands haye been getting deeper, while the amplitude 
of their diurnal yariation is diminishing, consistent with addi¬ 
tional CH 4 absorption at high northern latitudes rotating into 
yiew as the sub-Earth latitude moyes north. They conclude that 
the case for a yolatile transport contribution to the secular eyolu- 
tion looks strongest for CH 4 ice, despite it being the least yolatile 
of the three ices. 
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Table 7. Shifts (A) of the center of the absorption band of CH 4 . 


Longitude 

0.6190 pm 0.7296 pm 

0.7862 pm 

0.7993 pm 

0.8442 pm 

0.87 pm 

0.8897 pm 0.9019 pm 

49.7° 


-7.4 + 4.6 





-7.6+ 1.5 


103.7° 

- 

-15.6 + 2.7 

- 

- 

-9.3 + 1.7 

-13.9 + 2.2 

-13.9 + 1.8 

- 

141.3° 

- 

- 

- 

- 

- 

- 

- 

- 

197.9° 

- 

-12.1 + 1.4 

- 

-14.8+2.5 

-6.9 + 2.7 

-10.8 + 6.3 

-15.7 + 3.4 

- 

245.2° 

- 

-5.7+ 1.8 

- 

- 

- 

-5.5 + 2.6 

- 

- 

358.9° 

- 

-4.7 + 2.7 

-7.1 +2.0 

-4.0 + 3.0 

- 

-3.1 +2.6 

- 8 . 8 + 2.3 

- 


4.1. Longitudinal variations 

We note from Fig. |2] that spectra at different longitudes have 
different albedos as well as different overall shapes across the 
full spectral range, demonstrating the longitudinal heterogene¬ 
ity of the surface of Pluto. To show this in a more clear 
way, in Fig. |7] we plotted the variation with longitude of the 
spectral slope value obtained from our spectra, and compared 
it with a map of part of the surface of Pluto (mostly the 
northern hemisphere which is also the hemisphere currently 
observed from Earth) obtained by the New Horizon mission 
(www.nasa.gov/mission_pages/newhorizons). Notice that high¬ 
est (redder) value of the slope (at longitude ~ 100 °) is centered 
on the hemisphere that, according to the new maps of Pluto, 
contains the large dark region, whereas lower values correspond 
with a zone in average brighter (~240°). This is compatible with 
a higher abundance of solid organic material (redder in color 
and lower in albedo than the icy species observed in Pluto) in 
the darker region and a higher abundance of ices in the brighter 
ones. 

In Fig. 0 (bottom panel), we show the variation of the band 
depth of the five stronger CIT 4 absorption bands measured in our 
spectra. A sinusoidal trend is clearly visible with a minimum at 
a longitude between 140°-160° and a maximum at a longitude 
near 360°. This phase (360°) is also the one with the clearest de¬ 
tection of the 0.62-pm band. The location of this CH 4 maximum 
is in agreement with the trend found by G13 (maximum CH 4 
between 260° and 320°), but appears a bit shifted to the East. 
In particular, G13 finds that there are subtle differences in the 
longitude of the maximum of the absorption for different bands 
of CH 4 , with the weaker absorption bands having the minimum 
and maximum at longitudes more to the east than the stronger 
CH 4 absorption bands (see Eig. 4 of G13). They interpreted this 
difference between stronger and weaker bands as due to the pres¬ 
ence of different terrain types, with terrains with larger particles, 
or with greater CFI 4 enrichment, that produce the weaker bands 
more prevalent on the Charon-facing hemisphere. Eig. [ 8 ] shows 
how this trend continues for the bands in the visible that are 
even weaker than in the NIR, with the shallowest bands (0.80 
and 0.79 pm) having a minimum at ~180° and a maximum at 
~360°. The largest abundance of CFI 4 , or the largest particles, 
can be found in the hemisphere centered at ~ 0 °, and is not cou¬ 
pled either with the darkest (probably higher organic abundance) 
or with the brightest areas (probably higher volatile component). 

With respect to the degree of dilution of CH 4 in N 2 , G13 
found that all the CFI 4 bands in the infrared range show a blue 
shift (indicative of dilution), with a minimum blue shift cor¬ 
responding to the Charon-facing hemisphere, and a maximum 
blue shift corresponding to the anti-Charon hemisphere. That 
is exactly what Eig. |9] shows, with a maximum absolute value 
of the shifts in the anti-Charon hemisphere. The minimum of 
the shifts correspond with a region where the abundance of 



Longitude (degree) 


Fig. 8. Bottom Panel: variation of the band depth (in %) with longitude 
for the five stronger CH 4 absorption bands measured from the spectra 
of Pluto presented in this work. Top panel: variation of the depth of the 
absorption bands at 1.14, 1.72 and 2.32 pm from G13. The points in the 
gray shadow represent replicate values (at longitude ± 360°) to better 
visualize the variation of the band depth with the longitude. 


methane is higher (the Charon-facing hemisphere). This could 
indicate a higher degree of saturation of CH 4 in this h emisphere 
dTegler et alfcoiTUGrundv et alJl2013UTraftonl2015h . Our mea¬ 
surements in the visible show the same trend that the estimates 
from the NIR observations, apparently indicating that the degree 
of dilution is uniform over different path-lengths of the light. 


4.2. Secular variations 

In order to compare our observations with earlier sp ectral data, in 
Eig.flOl we show all our spe c tra, wi th spectra from iBuie & Einld 
(Il987 ) and iGrundv & Eint3 (Il996h overplotted on spectra ob¬ 
tained at similar longitudes (but more northern latitudes). Eor 
example, we compare our spectrum taken at a longitude of 49.7° 
(subEarth latitude 51°), with Grundy and Eink’s data for 1993 
(longitude 56°, subEarth latitude ~11)°, and our spectrum at 
358.9° longitude (subEarth latitude 51°) with 1983 spectrum 
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Longitude (degree) 


Fig. 7. Spectral slope (%/1000 A), calculated between 0.40 and 0.75 for spectra of Pluto normalized at 0.60 ^m vs longitudes 
(top panel), and map of part of the surface of Pluto result of combination of images taken by the Long Range Reconnaissance Imager 
(LORRI) on New Horizons, and combined with the lower-resolution color data from the instrument Ralph on the spacecraft (bottom panel) 
(www.nasa.gov/mission_pages/newhorizons). The scale of the horizontal axis is the same for the top and bottom panel. The red line indicate the 
subEarth latitude at which our data have been obtained (~51°). 



Fig. 9. Measured band shifts vs. Pluto’s longitude. Notice that at lon¬ 
gitudes between 100 ° and 200 ° the shifts are larger, thus the CH 4 :N 2 is 
smaller, which is compatible with a lower abundance of CH 4 suggested 
by the smaller band depths shown in Fig.[ 8 l 


from iBuie & FinkI (Il987h (longitude 329°, subEarth latitude ~ 


-10°). Differences in the new spectra compared to earlier data 
can be a result of the different latitude aspect, temporal changes, 
or a combination of both. From the comparision of the slopes in 
the visible (see Tabled, we can see that the spectra obtained in 
2014 are all redder than previous spectra. This could suggest a 
stronger presence of complex organics on Pluto’s surface with 
respect to previous epochs. 

G13 reported a difference in the secular evolution between 
stronger and weaker methane absorption bands observed in the 
near-infrared. They found that the stronger bands increase in 
strength with years, while the weaker bands in their sample do 
not. If this trend extends to the visible wavelength we would 
see a secular decrease of the depth of the bands. Indeed, from 
Figim we can note that the depths of the bands at 0.73 and 
0. 89 um measured in ou r spectra are lower than those reported 
in iGrundv & FinkI (Il996h . as we can expected for weak bands. 

G13 also found that the amplitude of the pattern of the lon¬ 
gitudinal variation decreases with time. This is consistent with 
the change over the years of the sub-Earth latitude of Pluto, with 
an increasing contribution to the observed spectrum of northern 
polar ices over the contributions from equatorial and southern 
latitudes, which show a much greater albedo contrast with longi¬ 
tude. In Fig. [TT] (top and middle panel) we plotted the measured 
band depth of the weak bands at 0.73 and 0.89 ttm, alon g with 
the depth of the same bands from iGrundv & FinM (Il996h . and a 
sinusoidal fit for each data set. In the bottom panel we show, for 
comparison purposes, the secular variation of the strong band at 
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Fig. 10. Secular variation on Pluto reflectance. Pluto-Charon spectra 
from this work shown with data at s imil ar spectral resolution an d signal 
precision from iBuie & FinkI dlQSTll an d I Grundy & FinkI lll996h at sim¬ 
ilar latitudes. All the spectra are normalized at 0.60 /rm and shifted by 
0.5 in relative reflectance for clarity. 


1.72 /mi from G13. The depths and the amplitude of the varia¬ 
tions measured by us are slightly lower than those of 1980 and 
1994, following the same trend of the near-infrared bands de¬ 
scribed in G13. 

In summary, the analysis of the secular variation of the slope 
in the visible, shows a reddening of Pluto’s surface from 1983 
to 2014, suggesting a higher degree of processing of the surface. 
On the other hand, the analysis of the secular variation of the 
methane ice bands in the visible is consistent with the results 
from the NIR, and indicates a decrease of the contrasts in band 
depth along Pluto’s surface. The clear correspondence of these 
changes with the pattern seen in the albedo maps of Pluto ob¬ 
tained with HST and New Horizons supports the contention that 
they are probably dominated by the sub-Earth latitudinal change 
rather than by seasonal changes. 

4.3. Comparison with Makemake and Eris 

To compare the spectral properties of Pluto, Makemake and Eris, 
in Sec.|3]we calculate for both the same parameters computed for 
Pluto’s spectra: slope and band depths. 


Fig. 11. Secular trend of the methane absorption bands: variation of the 
absorption band depth for the band at 0.73 pm (top panel) and 0.89 
pm (middle panel) compared with the secular variation for the band at 
1.72 pm (bottom panel) from G13. Lines in green represent a sinusoidal 
fit of band depths fro m this work; line s in re d represent a sinusoidal 
fit of band depths from iGrundv & PinB ( Il996lt . The points in the gray 
shadow represent replicate values (at longitude 360°) to better visualize 
the variation of the band depth with the longitude. 


The spectral slope of Pluto calculated between 0.57 and 0.70 
pm is identica l within the errors to that of Makemake from 
iLicandro et al.l (l2006bl) and three times larger than the spectral 
slope of Eris from Alvarez-Candal et al.l (1201 ih (see Table |3] and 
Fig- S- This suggests that the surface of Makemake contains a 
similar amount of solid complex organics (the coloring agent) as 
Pluto (averaging over the entire surface), while Eris’ surface has 
a smaller amount of organics th an both Pluto and M akemake, 
as previously suggested (e.g. bv iLorenzi et al.l[2015L and refer¬ 
ences therein). This shallower spectral slope is consistent with 
refreshing of the surface of Eris with volatile ices (N 2 , CO, CH 4 ) 
deposited by the collapse of its atmo sphere when th i s TNO was 
farther from the Sun, as sugg ested bv iLicandro et al.l (l2006ah and 
lAlvarez-Candal et al.l (1201 ll) . 


On the other hand, the methane bands are consistently deeper 
in the Makemake spectrum than in the Eris spectrum, and in 
the Pluto spectrum they are even weaker (see Table | 6 ] and Pig. 
Ell. This is indicative that methane is more abundant on the 
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Wavelength (^m) 


Fig. 12. Longitudinal variation inferred for MVICS observations. The 
average spectrum from the 2014 observations presented in this work 
is presented together with the filter transmission bands of the New 
Horiz ons Multispectral Visible Imaging Camera (MVIC) dReuter et alJ 
l2008h . The colored dots show the convolution of these transmissions 
with the average Pluto spectrum. The vertical bar shows the variation 
among the six reflectances (maximum and minimum value). The value 
at the NIR band is a lower limit as our spectra do not cover all the band¬ 
pass of that filter. 


surface of Maketnake than on that of Eris, and on the sur- 
face of Eris than on t hat of Pluto, as previously reported by 
iLicandro et al.l (l2006bl) . The two major known mechanisms that 
destroy methane ice on the surface of Pluto are transformation 
into red complex organics under energetic irradiation and subli¬ 
mation of methane into the atmosphere. As the overall colors of 
Pluto and Makemake are so similar, we infer that the more plau¬ 
sible explanation for the difference in the abundance of methane 
on Pluto and Makemake is the existence of a bounded atmo¬ 
sphere around Pluto. 

4.4. Other applications 

The spectra presented in this paper provide an independent cal¬ 
ibration of the MVIC (Multispectral Visible Imaging Camera) 
instrument aboard the New Horizons spacecraft, and allow the 
synthesis of a “green channel” that enable the preparation of 
false-color images of Pluto with the data from the spacecraft en¬ 
counter. Pig. [12] shows an average spectrum of Pluto from the 
present investigation with the transmission curves of the MVIC 
filters superimposed. Overplotted we show the photometric ap¬ 
pearance of this spectrum as seen through MVIC filters. The ver¬ 
tical lines indicate the maximum and minimum value of these 
“convolved photometric” points from our set of spectra and give 
an indication of how the variation showed by our data would be 
detected by MVIC. 

These new spectra provide the desired color information for 
modeling efforts that will use optical constants for non-ice mate¬ 
rials presently being synthesized in the laboratory by the UV and 
electron-irradiation of a Pluto mixture of ices known from near- 
IR spectroscopy (N 2 :CH 4 ;CO, 100; 1;1). The chemical analysis 
of the refractory residu e from a UV-irradiated Pluto ice mixture 
has been published by iMaterese et al.l (l2014l) . and results from 


the electron-irradiated ices are presented in Materese et al. (sub¬ 
mitted). 

5. Conclusions 

Low-resolution visible spectra of Pluto obtained at six different 
planetocentric longitudes between May 13 and July 17, 2014 
were obtained using the ACAM imager/spectrograph of the 4.2 
m William Herschel Telescope (Observatorio del Roque de los 
Muchachos, La Palma, Spain). These data cover a whole rotation 
of Pluto with a step of ~50°, and a gap of ~300°. 

The new spectra clearly show methane ice absorption bands 
at 0.73, 0.78 - 0.80 pm, and 0. 83 - 0.91 ttm also see n in previ¬ 
ously reported spectra of Pluto jGrundv & Einkl 19961) . We report 
the detection of the spectra of the methane-ice band equivalent 
of the 0.62-pm CH 4 gas absorption band in Pluto spectrum. This 
band was p reviously reported in the spectra of TNOs Eris an d 
Makemake (lAlvarez-Candal et al.ll201 ll ILicandro et al.ll2006bl) . 
The measurement of the depth of the absorption band of methane 
ice at 0.62 pm in the new six spectra of Pluto allowed us to es¬ 
timate the Lambert coefficient at this wavelength at T = 40° K 
(not measured in the laboratory yet) as 2.88x10^^ + 1.23x10“^ 
cm“^. Lrom the Makemake spectrum, we obtained a similar 
value within errors, 2.23x10“^ + 0.85x10'^ cm“'. Lrom the 
spectrum of Eris, we obtained a Lambert coefficient of 2.97x10“^ 
+ 1.46x10“^ cm“^ in this case for a T = 30° K. 

The spectral properties of Pluto, compared to those of Make¬ 
make and Eris, suggest that Pluto’s surface holds a similar 
amount of solid complex organics (the coloring agent) as Make¬ 
make, while Eris’ surface has a lesser a mount of organics than 


them both, as previously su 

igested (e.g. 

Licandro et alJl2006bl 

lAlvarez-Candal et al.ll201 U 

Lorenzi et all 2015). On the other 


hand, methane is more abundant on the surface of Makemake 
than on that of Eris, and on th e surface of Er i s than o n that of 
Pluto as previously reported by ILicandro et al.l (l2006bl) . 


The longitudinal variation of the slope in the visible is in 
phase with the distribution of the dark and bright material as 
seen in the albedo maps of Pluto from HST and New Horizons, 
with the reddest slopes in the hemisphere centered at 90° and the 
less red material in the hemisphere centered at 270°. 

The longitudinal and secular variations of the methane bands 
observed in the visible wavelengths are consistent with the be¬ 
havior of the weak absorption bands in the near-infrared region 
reported by G13. The maximum absorption due to methane ice 
is located at the phase centered at 360° (Charon-facing hemi¬ 
sphere) and the minimum at ~180° (anti-Charon-facing hemi¬ 
sphere). The curve is shifted from the curves of the stronger 
methane ice bands in the NIR. The depths and the amplitude 
of the diurnal variations of the methane ice in the spectra pre¬ 
sented here_(54^J9}are_slightly lower than those of 1980 and 
1994 dGrundv & Einkl [199^ . following the same trend of the 
near-infrared bands described in G13. This trend indicates that 
our observations capture regions with a lower contrast and is in 
agreement with the contrast in the maps of Pluto with more dark 
and bright material lying in the equatorial latitudes. 

All of the detected CH 4 bands are blue shifted, indicating a 
certain degree of dilution of the CH 4 in other ice, probably N 2 . 
We observe a longitudinal variation with a maximum shift lo¬ 
cated at longitudes between 100 ° - 200 ° and a minimum shift 
around 360° longitude. The fact that the shift is minimum in the 
regions where the abundance of methane is higher could be in¬ 
dicative of a dilution of CH 4 :N 2 more saturated in CH 4 . 
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